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Abstract. Stochastic dynamics of a vortex filament obeying local induced approx- 
imation equation plus random agitation is investigated by analytical and numerical 
methods. The character of a stirring force is supposed to be a white noise with spatial 
correlator concentrated at large distances comparable with size of the loop. Dependence 
of the spectral function ^s™s^ of the vortex line on both the one-dimensional wave 
vector « and intensity of the external force correlator (CkC«) was studied. Here s" is 
the Fourier transform of the line element position s a (£,t). It is shown that under the 
influence of an external random force a vortex ring becomes a small tangle whose mean 
size depends on external force intensity. The theoretical predictions and the numerical 
results are in reasonable agreement. 



1 Introduction 

In the previous paper Q we discuss how the large-scale perturbations can de- 
stroy the thermal equilibrium state in the space of vortex loop configurations. 
In this paper we elaborate that idea and present results of the both analytical 
and numerical investigations on stochastic dynamics of a vortex filament in Hell 
undergoing an action of the large scale random displacements. Moreover we con- 
sider the case when the smooth dissipation connected with normal component 
is small (that correspond to the case of very small temperature ) and the only 
strong dissipative mechanisms appear at very small scales comparable with the 
core radius of vortex. Thus detailed (at which scale) balance between the pump- 
ing and dissipation required for themal equilibrium is violated and, as it has 
been discussed in essentially nonequilibrium picture state must develop. 

From a formal point of view that problem is significantly more involved, 
therefore we restrict ourselves to a consideration of the local induction approx- 
imation in the equation of motion and omit processes of reconnection. 
This statement of problem is, of course, far from real superfluid turbulence in 
He II. A value of that work is that it enables us to understand mechanisms of 
entanglement of vortex filament and of appearing the strongly nonequilibrium 
state. We would remind that idea of the vortex tangle had been launched by 
Feynman more than 40 years ago Q and only about 10-15 years ago Schwarz 
demonstrated and confirmed that idea in his famous numerical simulations ||. 
To our knowledge a similar success in analytic study is absent. 

A structure of the paper is following. In the first part of this paper we de- 
velop nonequilibrium. diagram technique analogous to the one elaborated by 
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Wyld 0] for classical turbulence. Using further method of direct interaction ap- 
proximation we derive a set of Dyson equation for the pair correlators and for 
the Green functions. Assuming that region of stirring force and dissipation are 
widely separated in space of scales we seek for a scale invariant solution in the 
so called inertial interval. We also present results of direct numerical simulation 
of the vortex tangle dynamics. Numerical results confirm the ones obtained in 
the analytical investigations, however some discrepancies remained. 



2 Analytical Investigation 

In the local induction approximation (LIA) the equation of motion of quantized 
vortex filament in Hell reads 

a t 

Here s(£, t) is a point of the filament labeled by the Lagrangian variable £, < £ < 
which coincides here with the arclength. The quantity 5 stands for dissipation, 
which is small for usual scales and large for marginally small scales comparable 
with the core size ro. External Langevin force t) is supposed to be Gaussian 
with correlator 



<C a (6,*i)C /5 (£2,*a)) = - 6) - *a) ^ <*,p = 1,2,3, (2) 

where f a ((i-^) is changing on the large scale of order of the line length (~ 2ir). 

The quantity j3 is /3 = log ^ , with circulation K and cutting parameters 
R (external size, e.g. averaged radius of curvature) and r . For our numerical 

calculations we have chosen k = 10~ 3 cm 2 /s. This value corresponds to the 
case of superfluid helium. 

In Fourier space equation (Q) has the form 

- ^< = / r^lX-slSiq - qi - q 2 )d qi dq 2 + S q + Q ■ (3) 

Here is the spacial and temporal Fourier component of s Q (£,£), defined as 
follows: 

s g = J J sa (tt) e lM ~ K?) dtdi . (4) 
The vertex r"^ K2 responsible for nonlinear interaction has the form 

r^ R2 = ^Le^Kate-Ki), (5) 

where e"* 37 it the antisymmetric unit tensor. One can show the vertex r"^ K2 to 
satisfy the so called Jacoby identities 
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[K n r^ K2 + ,qrgg K1 + r& 3K ] s( K + Kl + « 2 ) = o, n = 2,4. (6) 

This relations express are tightly connected with the laws of conservation of total 
length L and curvature K 

L = I s 's ' d£ = const, K = I s "s " d(, = const . (7) 
Jo Jo 

Conservation of these quantities is readely derived from either of relations (JlJ) , - 
It is understood that conservation law is held in absence of the both dissipation 
and stirring force. 

One of the regular approaches to describe random fields is based on the Wyld 
diagram technique [Q, originally developed to study hydrodynamic turbulence. 
Following this technique we introduce for the description of random processes the 
following averages: the spectral density tensor (or correlator, or simply spectrum) 
S" & and the Green tensor G q P (or simply Green function) which are defined by 

Sf8{q + q x ) = (sX) . (8) 
/ Ss a \ 

Analysis of diagrams shows that due to the antisymmetry of tensor e Q ^ 7 con- 
tained in the expression for the vertex r"^ K2 , both S q P and G^ are propor- 
tional to 5 a p, i.e. Sg 13 = S q and G"P = G q . Details of that technique are 
described in ||. 

The renormalized quantities S q and G q (taking into account interactions) 
satisfy a Dyson set of diagram equations: 

Get opa , opa y^a/^ia ( ~\ C\\ 

q — + ^^"j i \ w ) 

S a q = G a q (F q a + *J) G q * . (11) 

Here °G q is the "bare" Green function which is equal to (oj — 5 re ) -1 . The mass 
operators <P q and S q can be written in form of diagram series: These series 
frequently used in nonequilibrium processes have a standard form, explicit form 
of them is given in H , ■ 

3 Conservation Laws and Pair Correlators 

Dyson equations have shapes indicating a cumbersome handling, therefore they 
can be studied for some special cases. One of them is considered in the present 
paper. It is connected with conservation laws for the total length and the curva- 
ture expressed by (0). Let us consider conservation of total curvature (for total 
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length there is the same consideration). In Fourier space the conservation laws 
for total curvature k can be expressed in the following form 

3K dP K 

where K K = s "s "e~ IK ^G?£ is the curvature density and P^ is the flux 

of this quantity in Fourier space (or, equally, in space of scales) . The right-hand 
side of equation describes creation of additional curvature (with rate I+(k)) 
due to external force and annihilation of it due to dissipative mechanism (with 
rate —I^(k)). In the equilibrium case the flux P^ is absent and source and sink 
terms must compensate each other locally for each k, i.e. I+(k) = I1(k). In 
the case under consideration when source and sink terms are widely separated 
in K-space that condition is obviously violated. Therefore a flux of curvature P^ 
in Fourier space appears. In region of wave numbers k remote from both region 
of the pumping k+ and of the sink k+ <C k <C k+ , the so called inertial 
interval, derivative dP^ /dn = , so P^ is constant equal ,say, P K . Resuming 



we conclude that the problem reduces to study the set of Dyson equation (10)- 
( pT| ) in inertial interval under condition of constant flux of the curvature. In this 
case Sg and G™ are expected to be independent on the concrete type of both 
the source and the sink but to be dependent on value of P K . Furthermore, the 
vertices r"^ K are homogeneous functions of its arguments. This property, as 
well as the condition k + <C k_ by virtue of which one can put k+ — and 
k_ = oo, leads to the assumption that the problem is the scale invariant, i.e. it 
has no characteristic scale for k. This suggests a power-law form of S 1 ™ and G™ 
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Sa = —ir/ — 1 . Go = —9 — ■ (13) 

q K r +P lK r l q K r L«fJ V ^ 

Here both / and g are dimensionless functions of their arguments. We aim now 
to find the scaling indices r and p. 

The first relation between indices r and p can be found from an analysis of 
diagram series, claiming all terms to have the same powers of argument k. This 
leads to the first scaling condition 



2r + p = 7 . (14) 

Another relation between r and p can be obtained from the Dyson equations 
([l0|), ( p"l| ) which can be rewritten in the form (see e.g. ||) 

J duIm{S«£2 - = 0. (15) 

This relation plays the role of kinetic equations for systems with a weak inter- 
action . It has been obtained with help of the expression for the Green function 
Gq = (u> — Sg)^ 1 ; the external force correlator F° disappears in the iner- 
tial interval. To find a relation of interest between r and p we rewrite relation 
(p"5|), disclosing expressions for mass operators 17™ and restricting ourselves 
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to first order terms in diagram series. That procedure called direct interaction 
approximation is frequently used in classical turbulence (see e.g. 0). After some 
calculation we arrive at the following relation (see also S): 



Im j dwdL0 1 dLJ2dnidK 2 5{q + q\ + q2) x r"^ 7 ^ 

X \ J KKiK 2 Lz <Z D qi°q 2 ' 1 k 2 kki Lt 92 D q D qi ' 1 «1 K 2 K Lz 9i °q 2 °q ) ~ V.[1V) 

To move further we perform conformal transformations in the second and third 
term within the braces in integrand, known as Zakharov transformations (see 
e.g. Zakharov [Q and Kuznetsov and L'vov ||). For example for the second 
term these transformations have the form 



- K "(k/k ) , Kl — K '(k/k ) , «2 = k(k/k ) , (17) 

u "(k/k "Y , ojx = u '(k/k ") r , uj 2 = u>{k/k "Y ■ (18) 



The third term is transformed in similar manner. As a result the integrand in 
16h becomes 



paj3-y fia r*/3 07 J 7-10/^7 ■ 
q qi q 2 ' KK -^ Kn 



K 2 K^i ' 



K 
Kl 



pf3*ya 



(19) 



where 



x = 7 - r - 2p . (20) 

Due to Jacoby identities (^) the integrand vanishes when x = — 2 for con- 
servation of total length and x = — 4 for conservation of total curvature. 
Substiting these values into ( |20| ) and solving equations (|l4|), ( |2C| ) we obtaine a 
set of couples of indices r, p corresponding to noncquilibrium. states with fluxes 
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of the length (r = 5/3 , p =11/3) and of the curvature (r = 1, p = 5). One 
time correlators can be found then integrating over frequences u> 



duj- 



1 



K r+P 



k " for length 
k~ 5 for curvature 



(21) 



So we have got solutions for the correlators S" which correspond to different 
conservation laws in (^) . Since there are no sources and sinks acting in the inter- 
mediate range these solutions guarantee that the according fluxes are constant. 
Depending on the way of agitation of the system one can get the real spectrum 
as some mixture of the obtained solutions in which the fluxes of length and cur- 
vature are present simultaneously. A similar situation for wave systems has been 
discussed earlier (see e.g. || and bibliography therein) and it is called multi-flux 
solution. 

Having in mind to compare our result with the both numerical and experi- 
mental investigation we have to take into account a presence of ^-correlated (in 
£-space) random force considered in our previous paper In the local induc- 
tion approximation the energy -ff{s} of line is proportional of its length and in 
parametrization when £ is arclength. can be expressed as 



47r r J 
r 



(22) 



It is easy to see that the equilibrium distribution described in [fi] leads in that 
case to correlator S" oc 1/k 2 . The final solution is a mix of equilibrium solution 
and of the ones expressed by relation (21 ). Because of nonlinearity it, in general, 
is not a simple superposition except of the cases when one of stirring action 
prevails and the other can be considered as small deviations. For instance if a 
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large-scale random stirring is small in comparison with (5-correlated (in £-space) 
action we have 

S " = ^ + ^ + ^- (23) 
The second and third terms in the right-hand side of ^ are small. The constants 
A, B, C entering are connected with both intensity of random stirring and its 
structure. The further specification requires some additional analysis. 



4 Some Numerical Results 

In this section we present some preliminary results on a direct numerical simula- 
tions of a vortex ring evolution under action of a random stirring displacements. 
The large scale character of noise was guaranteed by calculating it from a Fourier 
series taking into account only the first few harmonics. Besides some (uncon- 
trolled) white noise due to numerical procedure has been excited. Fig. [l] shows 
the projection of the line in the x, y - plane (where the ring was placed initially) 
for several times. As predicted, an consequent arising of higher harmonics takes 
place leading eventually to an entanglement of the initially smooth vortex loop. 

Another numerical results is shown in Fig. || where logarithm of quantity 
S% x averaged over several realizations is depicted as a functions of log k. The 
average slopes the graphs depend on intensity large-scale stirring force. In several 
realizations the slope lies between —2.5 and —3.5, which agrees with theoretical 
prediction ^3|. 

This work was partly supported under grant N 99-02-16942 from Russian 
Foundation of Basic research. 
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